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Microscopical model ol a doped lulleride electronic subsystem taking into account the triple orbital 
degeneracy of energy states is considered within the configurational-operator approach. Using the 
Green function method the energy spectrum of the model at integer band filling n — 1 is calculated, 
which case corresponds to AC@o compounds. Possible correlation-driven metal-insulator transition 
in the model is discussed. 

I. INTRODUCTION 

Electrical, optical and mechanical properties of fullerenes^ in condensed state demonstrate considerable physical 
content of phenomena which take place in fullerenes and show that the use of such materials in electronics has 
significant perspectives. Fullerene crystals and films are semiconductors with an energy gap of 1.2 — 1.9eV— & and have 
photoconductivity under visible light irradiation. Fullerene crystals have comparatively small binding energy and at 
room temperature the phase transition connected with orientational disordering of fullerene molecules take place in 
such crystals 5 . 

Addition of radicals containing platinum group metals 6 to fullerenes Ceo allows to obtain ferromagnetic material 
based on fullerene. In polycrystal Ceo doped by alkali metal superconductivity at temperature lower then 33A" is 
observed^. Large binding energy is typical for metallocarbohedrenes M§Ci2, where M = Ti, V, Hg, Zr. For example, 
in Ti 8 Ci2 molecule binding energy per atom is 6.1eV— (for Ceo molecule this energy is 7.4 — 7.6eV— ). 

Fullerenes in solid state (fullerites) are the molecular crystals, where interaction between atoms in Cqo molecule is 
much larger then interaction between nearest molecules. In tightly packed structure each fullerene molecule has 12 
nearest neighbors. Depending on peculiarities of molecular interaction, face-centered cubic lattice or hexagonal lattice 
is realized^. Phase transition in Cqq crystal occurs at the temperature of 257A and this is the first order transition. 
At high temperatures molecules can freely rotate whereas at low temperatures rotation is stopped and anisotropy 
of neighbor molecule Ceo interaction becomes important. This leads to small sharp change of distance between the 
nearest molecules. According to results of X-ray structure analysis^ lattice constant changes from 1.4154 ±0.0003nm 
to 1.4111 ± 0.0003nm (that is by 0.43 ± 0.06 percent). 

At low temperature, when Cgo -molecules are oriented in space, crystal lattice symmetry does not coincide with 
the symmetry of single molecule Cqq (icosahedral symmetry Y) . In a unit cell of fullerite crystal lattice there are four 
Cgo-niolecules. These molecules form tetrahedron in which orientations of all molecules are the same. Tetrahedra, in 
their turn, form simple cubic lattice. 

Fullerites are semiconductors with energy gap of 1.5 — 1.95eU 3 . Electrical resistivity of polycrystals Cecr^ mono- 
tonically changes with changing temperature and energy gap has monotonic dependence on the pressure value: an 
increase of energy gap under the pressure, higher than 2 x 10 5 atm indicates the absence of metal-insulator transition 
at p ~ 10 6 atm. In the temperature region 150 — 400A the relaxation time is temperature-independent what indicates 
that the carriers are localized and hopping mechanism of recombination, which includes tunneling of electrons between 
localized states, is realized. 

It has been shown in 1991— that doping of solid fullerenes Ceo by small quantity of alkaline metal leads to formation 
of material with metallic type of conductivity and this material becomes superconducting at low temperatures (T c from 
2. 5 A for Na2KC6o to 33 A for RbCs2Cgo)- At changes of temperature, concentration of alkaline metal, parameters 
and structure of lattice various phases of these compounds have been realized. In particular, at various filling n (n 
may change from to 6) of lowest unoccupied molecular orbital (LUMO) the metallic, insulating or superconducting 
phases have been realized. Superconductivity in doped fullerenes K^Cgo has been studied theoretically in paper— and 
strong electron correlations have been shown to play a crucial role in superconducting state stabilization. Recently, 
strong electron correlation were also proverb to be responsible for superconductivity of planar carbon systems of 
graphene type. 

Let us consider the electronic structure of Ceo in detail. In single-particle approximation, neglecting electron 
correlations, the following spectrum has been calculated^: 50 of 60 p z electrons of a neutral molecule fill all orbitals 
up to L = 4. The lowest L = 0,1,2 orbitals correspond to icosahedral states a g ,ti u ,h g . All states with greater 
L values undergo the icosahedral-field splitting. There are 10 electrons in partially filled L = 5 state. Icosahedral 
splitting (L = 5 — > h u + ii„ + t^u) of this 11-fold degenerate orbital leads to the electronic configuration shown below. 
Microscopic calculations and experimental data show that the completely filled highest occupied molecular orbital is 
of h u symmetry, and LUMO (3-fold degenerate) has t\ u symmetry. At such conditions HOMO-LUMO gap appears 
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due to icosahedral perturbation in the shell with L = 5; energy gap found experimentally is about leV for molecules 
in vacuum. A t 2g (LUMO+l)-state, originated from L — 6 shell, is found approximately leV above the t\ u LUMO. 

Electron-electron correlations in Cgo are described by two main parameters: intra-molecular Coulomb repulsion 
U and Hund's coupling Jh- In fullerenes the competition between intra-site Coulomb interaction (Hubbard U) 
and derealization processes, connected with translational motion of electrons (which determines the bandwidth), 
causes the realization of insulator or metallic state**. Majority of the experimental data and theoretical calculations 
indicate that all materials with ions C^ ™ at integer n are Mott-Hubbard insulators as U is quite large for all doped 
compounds A^Cgo. Fullerides A x Cqq doped with alkali metals A attract much attention of researchers due to unusual 
metal-insulator transition in these compounds. Only A3C60 is metallic and other phases ACgo, A2C60 and A4C60 
are insulator 15 . This experimental fact contradicts to the results of band structure calculations (seei£ for example) 
which predict purely metallic behavior. It has been noted in paper 17 , that for explanation of metallic behavior of 
Mott-Hubbard system [x = 3 corresponds to the half-filling of the conduction band) one has to take into account a 
degeneracy of energy band. On the base of Gutzwiller variational approach the metal-insulator transition has been 
proven 1 ! to exist for all integer band fillings. It is shown that the critical value of Coulomb interaction parameter 
depends essentially on the band filling and degeneracy (in case of half filling ~ 2,8 for double de generacy, 
~ 3,9 for triple degeneracy). The present study is devoted to investigation of Mott-Hubbard localization in 
electronic subsystem of fullerides with strong electron correlations within the model taking into account the orbital 
degeneracy of energy levels, strong Coulomb interaction and correlated hopping of electrons. 

II. THE HAMILTONIAN OF DOPED FULLERIDE ELECTRONIC SUBSYSTEM 

Within the second quantization formalism the Hamiltonian of interacting electron systems can be written^ as 

i\a ijA<T ijkl afiyS erer' 



where the first sum with matrix element 



Uj = J d 3 rcj)* Xi (r - Ri) x 



2 



A + V ton {r) 
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(r - R,-) (2) 



describes translational motion (hopping) of electrons in the crystal field V lon (r) and the second sum is the general 
expression for pair electron interactions described by matrix elements 

JtfkT = J J^ T - R ^( r - Rtf) x rr=77i ~ R 0M r ~ Kk)drdr>. (3) 

In the above formulae a^ X(J , are operators of spin-cr electron creation and annihilation in orbital state A on lattice 
site i, respectively, indices a, /3, 7, 5, X denote orbital states, <p\i is wave-function in Wannier (site) representation 
other notation are standard. Hamiltonian ([1]) is essentially non-diagonal and hard to treat mathematically. The 
problem can be greatly simplified by neglecting the matrix elements of interaction of the third and further orders 
of magnitude and restrict oneself to consideration of a single orbital per site. In this way, Hamiltonian of Hubbard 
model and many other backbone models of strongly correlated electrons theory were derived. However, it has been 
shown that these models lack the possibility of description of electron-hole asymmetry, observed in real correlated 
electron systems. To maintain such possibility we are to consider the energy levels structure and estimate interaction 
parameters prior to make simplifications. Following paper o 19 i 20 we derive the Hamiltonial which takes into account 
the correlated hopping of electrons (the site-occupation dependence of hopping parameters results from taking into 
account the interactions with second order of magnitude matrix elements) and variety of intra-cite interactions caused 
by triple orbital degeneracy of LUMO in doped fullerites. Interaction integral of zeroth-order magnitude is on-site 
Coulomb correlation (characterized by Hubbard parameter U): 



U = J J |0*(r-R 4 )| 2 ^— —\^{v' -Ri)\Hrdr' , (4) 
In orbitally degenerate system, the on-site (Hund's rule) exchange integral 

Jh= f I <«(r ~ Rt)4> x > (r - Ri) I — —rr x 0* v (r - R<)0 A (r' - R z )drdr , (5) 
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is of principal importance, too. Parameter U value for fullerenes have been estimated within different methods. Use 
of local density approximation (LDA) gives 3.0 e V 21 ' 22 . Experimental estimation of electron repulsion energy^ 3 - gives 
C/~2.7eV. 

It's worth to note, that in solid state molecules are placed close enough to provide substantial screening of interaction. 
Calculation with screening effect took into account give U 2.7 21 i 22 . Combining Auger spectroscopy and photoemission 
spectroscopy lead to value 1.4-1.6 e V 24 ' 25 for U. We also note that energy cost of electron configurations with spins 
aligned in parallel is considerably less than for anti-parallel alignment. Orbitally degenerate levels are filled according 
to Hund's rule. Experimental methods 2 ^ for singlet-triplet splitting give 0.2 eV ±0.1 eV; and in work- 2 ^ has the values 
close to 0.05 eV. The relevant inter-site parameters are electron hopping integral and inter-site exchange coupling 
J(i\jX'jXi\'). 

The resulting Hamiltonian of doped fullcridc electronic subsystem reads as 



i\a iX iXrr iXX'a 



U' U' - J H 
~~2 n i^ni\ s H ^ 

iXa 

+ Y M n K + A<x a jA CT + Y 4 { a t\<T a 3>>° n i\ + h - c ) + Y 4 {a+^ajxan^g + h.c.) , (6) 

ijXa ijXa i'jXa 

where riix a = o-txa a ^<r , U' = U — 2Jh and hopping integrals (n) , , taking into account three types of correlated 
hopping of electrons 2 ^ are introduced. 

In a model of triply degenerate band, every site can be in one of 64 configurations (see fig. 1). To pass from electron 
operator to Hubbard operators X pl of site transition from state \l) to state \p) we use relations of type 

^+ _ TOO, 000 _|_ ^200,4-00 _|_ ^tt0,0f0 _|_ j^ti0,0t0 _|_ j^-fOt.OOt _|_ j^-tOJ.,004. _|_ j^t20,020 _|_ j^t02,002 

+ ^24.0,440 + X 2oum + _x"rtt,ott + x niMi + j^2to.ito + X 2nt,m + x ni,on + x m,oit 

+ ^244444 _|_ j^24t,44t + j^2t4-4t4- + ^2tt414 + ^220420 + x f2t,02t + ^24,024 + ^202402 

_|_ j^-tt2,0f2 , x tl2M? _i_ j^t22,022 , j^224424 , j^242 442 , j^22f42f , J^2f24t2 , j^-222,4-22 
-+ _ j^400,000 _ ^200400 _|_ j^4t0,0t0 _|_ j^440,040 _|_ j^40t,00t _|_ j^404,004 _|_ j^420,020 _|_ ^402,002 

_ ^2f04to _ j£-20t,f0t + X m,on + 444,044 _ ^240440 _ ^204404 + x iUfin + _^44t,04t 

_ X 2n-m _ X 2n,m _ ^244444 _ ^244444 _ ^220420 + ^42t,o2t + ^424,024 _ ^202402 

_|_ j^4T2,0t2 _|_ j^442,042 _|_ j^422,022 _ jj£-22t42f _ J^2f2,tt2 _ j^224424 _ j^242442 _ 222422 fy\ 

which ensure the fulfilment of anticommutation relations {Xf ; Xj 1 } — 5ij{5ikXf f + S p tX^ L ), and normalizing condition 
Xf = 1, for number operators Xf — X? X? of \p >-state on site i. Such type of electronic operators representation 

i 

is typical for models of strongly-correlated electron systems as superconducting cuprates^ 9 -, manganites 3 - ., cobaltites 3 ^, 
optical lattices 3 - 2 -! 3 * 3 . Using the root vector notations introduced in paper-24 allows to obtain much more compact form 
of Hamiltonian in configurational representation. However, in our case number of subbands is relatively small and we 
use bulky but simple notations which make the projection procedure used below more transparent. 

In the configurational representation the model Hamiltonian takes the form H = Hq + T. Here Hq sums the "atomic 
limit" terms and the translational part may decomposed as T = ^ T nm , where n,m serve for numbering "atomic" 

states. Terms T nn of the Hamiltonian form the energy subbands and terms of T nm describe the hybridization of 
these subbands. Different hopping integrals correspond to transitions in (or between) the different subbands. The 
subbands of higher-energy processes appear to be narrower due to the correlated hopping of electrons. The relative 
positions and overlapping of the subbands depends on the relations between the energy parameters. At integer values 
of electron concentration (n = 1,2,3, 4, 5) in the system the metal-insulator transition is possible. 

In the partial case of band filling n = 1, strong Coulomb correlation and strong Hund's coupling (parameter U — 3Jh 
is much greater than the bandwidth, see estimations in paper o 21 ' 22 ) the states with three and more electrons on the 
same site are excluded. Then the influence of correlated hopping can be described by three different hopping integrals. 
The bare band hopping integral is renormalized to take into account the band narrowing caused by concentration 
dependent correlated hopping as Uj(n) = Uj(l — Tin). This hopping integral characterizes lower Hubbard subband. 
Parameter t± is usually neglected, but it is of principle important for a consistent description of correlation effects 
in narrow band systems (se o 19 ' 20 for a detailed discussion). The hopping integral for upper Hubbard subband is 
Uj(n) = Uj(n) +2iy- and fy(*n) = tij(n) +ty- describes a hybridization of lower and upper Hubbard subbands. In the 
following only the case n = 1 is considered so we omit the explicit notation of concentration dependence. Then the 
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FIG. 1: Possible site configurations in threefold degenerate model. The first symbol in the state notation correspond to a 
orbital, the second and the third - to /3 and 7 orbitals, correspondingly. 



Hamiltonian in X— operator representation^! has the form 
H =H Q 



(A) , TfW 

b + H h 



(8) 



A— a,/3,7 
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rj(a) V~V+ vcrOO.000 v OOO,erOO , 7 v crcrO.OcrO v Oa0.craO . 7 v aOa.OOa v OOa.crOa 
<<■, ZJ''-' • A , + % A i A , A , 

ija 

+ HjX i A j + Wj" A i A 3 J) 

Tj(a) Y^Z /^crOO.OOO v 0cr0,crcr0 . ^crcrO.OcrO v 000,cr00 , 

v ctOO,000 v OOa,nOa , v aOcr,OOa v 000,o-00\ 

it(P) ^OerO.OOO ^OOO.OerO . 7 v aaO,<rQO v aOQ.aaO , 7 v Oaa.OV)(J V C 

2J'</ A . . X j + t ij A i A , 



A , A .; _ '<.; A , A ., J 



zt(/3) Y^T I v OcrO,000 v 00ct.0ct(t . v 0erer,00o- ^OOO.OctO v OctO,000 ^ctOO.ctctO ^o-o-O.erOO v OOO.OctO\ 
ija 

''/, ZJ''-' A < / A j A , A , + 



H h =-2_J lJ {X l Xj + X t Xj +X Z Xj +A, Xj ). 



Green functions technique allows us to calculate the energy spectrum of the model which corresponds to the 
electronic subsystem of A^Ceo in the case of electron concentration n = 1. One can rewrite the single-particle Green 
function ((aiXal^x^)) on the basis of relation between electronic operators and Hubbard's X-operators: 

a pat = AT p 00( W° + X?°>™ + x oot,tot = X ooo, t oo + ^ (9) 

where the operator Y p describes the transition processes between doubly occupied Hund's state and single occupied 
state. The processes involving other type of doubly occupied states, empty states, states with three or more electrons 
is improbable due to energy scaling. 

In this way we obtain the following expression for the single electron Green function 

««wi4«t» = {(^ 000 ' t00 i4° ' 000 )) + +<<*r t00 i^>> + <<^r t00 >> + «W »• ( 10 ) 

Equation of motion for Green function ((^° 00 ^ 00 |^ no ' non )) has the form 

vooo I -rtoo 

(E + ^iix^ixl? ' 000 )) =v — * +(([^ 000 ' T00 ;E^ A) ]l A 'p° 0,000 )) 

A 

+ (([X° 00 'too.^ i/ (A) ]|x Too,ooo )) (n) 

A 

and equation of motion for Green function ((^plA'p, 00 ^ 00 }} - 

(E + M - U + 3J H )((Y p \X°;°^)) = (([Y p ; £ H™] l^ 000 )) + ((Y p , £ H^]\X^ 000 ))- 

A A 

To obtain closed system of equations for Green functions (( X° 00 ^ 00 \ X^, 00 ' 000 )) and ((Y^X^ ' 000 }} we use the 
projection procedure similar to the worker 



[* P 000 ' t00 ;E^ A) ] = E4^r' t0 °; (12) 



jW] _ \^,J> v 000,t00 
A i 

A i 

K>;E^ A) l - E4^; 
[^;E^ A) ] = E^ 000 ' t0 °. 
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As a result after Fourier transformation we obtain the Green function in the form: 

{{Xi ^ ))k ~ 2^ X (E - El (k))(E - E 2 (k)) > (13) 

where the quasi-particle energy spectrum 

6/V\ _i_ ?bf 



E lt2 (k) = - M+ IL^lE + ! W+ £ W T ^(f/ _ 3J ff - e *(k) + e»(k)y + 4 £ "(k)s"(k). (14) 

In the absence of orbital order the energy spectrum for f3 and 7 electrons is the same as for a electrons. 

The non-operator coefficients e (k) , e b (k) , e h (k) , e h (k) one can obtain by the anticommutation of Eq. (Ti"2l) with basis 
operators x® 00 '^ 00 and and following replacement of operators by c-numbers (see in this connection^) . 

4 = ^[ik(<* p 000 (^r + 4, 00 )) + +(K m ( X ° P ' Q + K™» + (Xf * 00 X^ 00 ) + (X oto )t oo x too,oto ) + 
_l_^oj.o,too^too,o4,o^ _|_ _|_^oot,too_ytoo,oot^ _|_ ^oo4,,tooj^too,oo4.^ {^((X^'^X 000 ^ ) + (X tot,ooo x ooo,tot^j 

£ h = -Lt k l((X°™ + x; 00 ) x (X°, t0 + X°, ot + A™ + #)) + (XOtO'tooxtot.ott) + (xn o - 000 x° 00 'ttO) _ 
_^oto : too^too,oto^ _ ^oot,too^too : oot^ _|_ ^tpt.ooo^-ooo.tot^ _ ^OOt.tOOj^ttO.Ott^j 

4 = -^[(xtto.ooo^ooo.no) + {x ^ 000 x a p oa ^)] + 

+^[(^p 0t0 + + X^ + Xt°t) x (x o;o + x m + x oo t + jftpt) + (x oto,too x too,oto ) + ^ott.tot^tot.ott) _ 



_ (x oto,too x tot,ott ) _ (x on,tot x too,oto ) + (x oo t , t oo x too,oot ) + (x on,no x tto,on ) 

+ (x 00 ^^ 00 x no ' on ) + (x° tt ' tt °X 1 ' 00 ' 00t )] 



c [((X° p ?° + Xl?°) x (XW + x tot + x o t o + x oo t)) + {x ott,tot x too,o t o ) _ (x o t o,too x too J o t o ) + 
+ (x tto,ooo x ooo,tto ) + +(x oot,roo x too,oot ) _ (xW' 000 X p 000 ^ t)(X p °tt.tto x Too,oot )]j 

where C x = (X°°° + A^ 00 ), C 2 = (X^° + X^ + X^° + if), it i s wor th to note that in the partial case of band 
filling n = 1 and strong Coulomb correlation we work with reduced Hilbert space of electronic states, so C\ + C 2 = 1. 

Let us denote the concentration of empty lattice sites by e, concentration of singly occupied sites with spin a 
electron in orbital state A by s\ a , Hund's doublons concentration by d a , Hubbard doublons by d 2 and non-Hund 
doublons by d. In a paramagnetic state s\ a — s, d a — d. For the case of strong Hund's coupling the high energy 
doublon configurations are excluded, d 2 — d = 0. We can utilize the completeness condition for the X-operator set 
to have constraint e + 6s + 6d = 1, which, at condition e = 6d, leads to the equation 

• = — e— ( 15 ) 

Finally in the paramagnetic case at n = 1 we obtain 

h 216d 2 -12d+l 72d 2 - 

e= 24rf+ i tk + 2im fk; (16) 

, Id - I2d 2 

£ =ik ^6^' (17) 

h 36d 2 i k , . 

£ ^ k l-6rf + 2(T-6d)' (18) 

, 24d + 1 - 216d 2 , x 

g =ffc 3(24^+1) > ^ 



7 



In this way, the energy spectrum depends on the concentration of doublons d (through the dependence of non- 
operator coefficients). The doublon concentration is determined by the condition 

6 j- 1 W Mk) I Be(k) ^ (20) 

where 

U - 3 J H + e b - s b \ 




^{U - ZJ H -e b + e b ) 2 + Ae h e h J ' 
A e (k). (21) 

Using the model rectangular density of states at zero temperature one obtains 

twj-v E-E x {e) 4wJ_ w E - E 2 {e) 

here 0(— E(e)) is Heaviside theta-function. Solving this equation numerically we obtain the doublon concentration 
as function of the model parameters. To study a metal-insulator transition (MIT)2& — we apply the gap criterion 

AE = E 2 (—w) - Ei(w) = 0. (23) 

In the point of MIT the polar states (holes and doublons) concentrations equals zero. Thus, for the non-operators 

2 ' fc — 3 

AE = U - 3 J H - w - w. (24) 



coefficients we have e b = ik, £ h = 0, e b — 4jr, e h = and for the energy gap we have the equation 



Here w = z\t\(l — n) and w — z\t\(l — ti)(1 — 2r) are the halfbandwidths of the lower and upper subbands, 
respectively, z is the number of nearest neighbours to a site, \t\ is the magnitude of bare nearest-neighbour hopping 

integral, t±,t = jr^-; are the correlated hopping parameters. From the equation (|24l) one obtains that the critical 
value of the intra-cite Coulomb interaction parameter equals the sum of quasiparticle subbands halfbandwiths. 

Analysis of the expression (l24l) allows explaining the differences of electrical characteristics (insulator or metallic 
state realisation) depending on the correlated hopping strength. 

The correlated hopping influence substantially on electrical characteristics of narrow band material with three-fold 
orbital degeneracy of the energy levels. Both the filling of the sites involved into the hopping processes (through the 
correlated hopping of the first type) and the neighbor sites (through the second type correlated hopping), can lead to 
appearance of the gap in energy spectrum and stabilization of the insulator state. The energy gap, however, opens 
at relatively large increase of correlated hopping parameters which can not be achieved in a compound by change of 
external conditions only. Such critical increase of parameters T\ and r can be realized at doping. A distinct picture is 
observed at the change of intra-site Coulomb interaction parameter. At increase of (U — 3Jh)/w over a critical value 
(dependent on the correlated hopping strength) the energy gap occurs and the metal-insulator transition takes place. 
The critical value for the partial case of the model when the quasiparticle subbands have the same widths (in absence 
of the correlated hopping), is (U — 3Jh)/w — 2 which corresponds to the result of work a 28 i 38 for non-degenerated 
Hubbard model. 



III. CONCLUSIONS 



Within the variant of triple orbitally degenerate model of the electronic subsystem of a doped fulleride compound 
considered above not only the on-site Coulomb correlations but also additional interactions of principal importance, 
namely the correlated hopping, can be introduced and analyzed. The use of Hubbard X-operators representation 
appears to be useful to exclude from consideration the parts of Hilbert space which are irrelevant at particular band 
filling. The ground state metal-insulator transition in the triply degenerate model of partially-filled doped fulleride 
band takes place at moderate values of the correlation parameter which in this case is a combination of on-site Coulomb 
repulsion energy, Hund's rule coupling and electron hopping parameters. The correlated hopping of electrons leads 
to further localization of current carriers. The influence of the correlated hopping is substantial and makes the 
estimation of the model parameters from the available spectroscopic data ambiguous. The problem can be resolved 
by the additional spectroscopic experiments with use of external pressure. In this case the reasonable estimates could 
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be obtained using the fact that in distinction from the on-site parameters, the correlated hopping parameters must 
be pressure-dependent. The metal- insulator transition described above can be realize d 39 ' 40 in the doped fulleride 
compound under the external pressure. 
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